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摘要 :植被 指数 是 研究 区 域 植被 变化 的 重要 手段 。 基 于 MODIS/Terra NDVI 遥感 影像 数据 ,对 陕西 榆林 市 2000 一 2016 年 8 月 份 
NDVI 进行 像 元 尺度 的 时 空 变化 及 影响 因素 分 析 。 结 果 表 明 ,榆林 有 96.44% 的 区 域 植被 指数 是 增加 的 3 增加 值 在 0 一 0.02/a 区 
域 面积 占 榆林 的 93.63% , 呈 显 著 增 加 趋势 的 面积 占 榆林 的 80.72% 。 复 直线 回归 分 析 显 示 ,气象 因素 对 植被 生长 ,演化 起 到 了 
促进 作用 ,人 类 活动 使 得 植被 指数 等 级 变 得 更 加 均匀 。 人 类 活动 抑制 植被 生长 演化 的 区 域 占 榆林 总 面积 的 45.04% , 主要 分 布 
在 榆林 市 最 北部 的 府 谷 县 榆林 市 南部 及 西部 的 大 部 分 区 域 ;人 类 活动 促使 植被 指数 增加 的 区 域 占 榆林 总 面积 的 S4.96% ,说 明 
有 一 半 以 上 的 区 域 , 人 类 活动 对 植被 生长 起 到 了 促进 作用 ,这 些 地 区 的 封 而 育林 \ 退 耕 还 林 < 退 牧 还 草 等 措施 的 实施 效果 较 好 。 
关键 词 :植被 指数 ;时 空 变化 ; 复 直 线 回 归 ; 像 元 尺度 ;榆林 
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Abstract: Vegetation index is .an important method to study regional vegetation changes. The present study is based on 
MODIS/Terra NDVI remote sensing data, and analyzes the spatial-temporal changes and impact factors of the NDVI in 
August from 2000 to 2016 at pixel scale in Yulin City of Shaanxi Province. The results showed that vegetation index 
increased' in approximately 96.44% of the land area in the Yulin city，and the area with increasing NDVI between 0 一 
0.02/a accounts for 93.63% of the land area of Yulin city. The area with significant increasing trend accounts for 80.72%. 
The complex linear regression analysis showed that meteorological factors play a positive role in the growth and evolution of 
Vegetation ，and human activities also make the vegetation index to become more uniform. The area where the human 
activities restrain the growth and evolution of the vegetation covers 45.04% of area of Yulin. lt is mainly distributed in Fugu 
County, which is located to the north of Yulin and most areas of the southern and western parts of Yulin. The area where 


human activities promote the increase of the vegetation index accounted for 54.96% of area of Yulin which indicated that in 


基金 项 目 :陕西 省 自然 科学 基础 研究 计划 (2016JQ5105 ) ;榆林 市 科技 计划 项 目 (2016-21) ;国家 自然 科学 基金 项 目 (51609270 ) ;教育 部 博士 学 科 
点 基金 (20136118120021) ;西安 理工 大 学 水 利水 电学 院 科研 项 目 (2016ZZKT-15) 
收 稿 日 期 :2017-04-21; 网 络 出 版 日 期 :2018-01-08 


* 通讯 作者 Corresponding author. E-mail: liudf@ xaut.edu.cn 


http :// www.ecologica.cn 


1 AVIAN 人 下 半日 二 | 
haX iV 合作 期 刊 


站 hh 
\ | 
Hl 


8 期 栾 金山 等: 近 17 年 陕西 榆林 植被 指数 的 时 空 变化 及 影响 因素 2781 


more than half of the region, human activities have played a positive role in the growth of vegetation. In these areas, the 
effects of forest conservation, returning farmland to forests, and returning pasturage to natural grassland are greater than 


those in other areas. 


Key Words: vegetation index; spatial-temporal change; complex liner regression; pixel scale; Yulin 


地 表 植 被 处 于 岩石 圈 最 活跃 的 一 层 ,是 岩石 圈 与 大 气 圈 、 土 壤 圈 、 生 物 圈 水 圈 长 期 相互 作用 的 产物 ,在 地 
球 生态 系统 平衡 .气候 变化 和 水 循环 中 起 着 协调 的 作用 '" 站 i。 同 时 ,植被 是 全 球 土地 覆盖 最 重要 .最 敏感 的 构 
成 要 素 之 一 ,可 在 较 短 时 间 内 反映 气候 变化 和 人 类 活动 的 影响 5 。 土 地 沙漠 化 敏感 区 域 属于 全 国防 风 固 沙 
生态 功能 区 ,对 国家 生态 安全 具有 重要 作用 ,关系 着 区 域 人 地 了 矛盾 的 解决 的 。 榆 林 位 于 毛乌素 沙漠 与 陕 北 黄 
土 高 原 交 界 地 带 , 属 于 农 牧 交错 区 ,生态 环境 脆弱 上 且 变 化 敏感 ,是 研究 沙漠 的 热点 和 重点 区 域 [551 : 

植被 指数 时 空 变化 影响 因素 的 研究 大 多 以 相关 分 析 为 基础 。Gao 等 "认为 在 中 国 东 部 ,气温 和 降水 同时 
影响 植被 变化 ,但 气温 起 主导 作用 ,同时 得 出 NDVI 对 气温 和 降水 的 变化 有 滞后 期 ;分 别 为 104 和 30d。 也 有 
些 学 者 分 析 后 发 现 NDVI 对 降水 及 温度 的 响应 具有 明显 的 空间 差异 与 滞后 效应 人。Piao 等 :59 应 用 相关 分 
析 方 法 分 析 了 草地 NDVI 气候 变化 响应 。 夏 照 华 '" 分 析 了 整个 中 国 最 近 宇 十 多 年 来 NDVI 的 时 空 分 布 规 律 ， 
以 及 NDVI 与 降水 、 温 度 之 间 的 关系 ,尝试 区 分 出 了 NDVI 影响 因子 中 的 气候 因子 与 人 类 活动 因子 。 
Fabricante 等 2 认为 NDVI 值 的 变化 与 前 几 个 月 的 水 总 量 相关 更 大 。 Ding 等 利 用 NDVI 得 出 青藏 高 原 北 
坡 的 生长 期 只 有 3 个 月 (7 一 9 月 ) ,生长 期 的 降水 与 NDVI 的 相关 性 极 好 。 李 登科 等 :外 利用 简单 相关 系数 分 
析 , 分 析 了 长 城 沿线 各 区 县 年 均 NDVI 与 降水 和 气温 的 相关 关系 杨 英 莲 '5 分 析 了 青海 省 不 同 草地 类 型 
NDVI 的 季节 变化 和 年 际 间 的 变化 特征 ,又 将 22 年 来 草地 平均 NDVI. 气 温和 降水 进行 了 趋势 分 析 和 相关 分 
析 。 李 丽 娜 "在 空间 中 选取 足够 多 的 点 对 陕西 省 生长 季 NDVI 与 温度 .降水 进行 相关 分 析 , 发 现 不 同 区 域 气 
候 对 植被 影响 有 差异 ,植被 对 气候 的 敏感 度 也 有 差异 。 朱 文彬 等 '”) 通过 逐 像 元 分 析 发 现 柴 达 木 盆地 植被 覆 
盖 的 影响 因素 主要 包括 降水 .地 表 地 下 水 文 条 件 海拔 高 度 和 人 类 活动 4 个 方面 。Zhao 等 ' 对 1982 一 2003 
年 新 疆 NDVI 进行 研究 发 现 ,区 域 NDVI 的 增加 与 降水 量 和 潜在 蒸 散 量 (ET) 的 增加 有 关 。 徐 浩 杰 等 :2 利用 
简单 相关 系数 分 析 ,发现 影响 祁 连 呈 植被 生长 的 主要 因子 是 气温 和 降水 ,局 部 地 区 人 类 活动 也 是 影响 植被 生 


长 的 关键 因子 。 张 清 雨 等 .站 对 内 蒙古 自治 区 的 NDVI 进行 研究 和 分 析 , 发现 整 个 自治 区 大 部 分 区 域 NDVI 与 
年 降水 量 呈 显著 相关 ,同时 还 发 现 近 :30 年 来 人 类 活动 对 植被 NDVI 的 影响 程度 正在 逐渐 增强 。 张 智 炸 等 1 
对 影响 大 豆 NDVI 的 气象 因素 进行 多 元 线性 回归 分 析 , 从 而 实现 对 NDVI 的 预测 。 孙 庆 龄 等 ' 引 应 用 逐 像 元 的 
气候 因子 与 生长 季 累 计 NDVI 的 相关 性 分 析 的 方法 ,分 析出 影响 三 江 源 地 区 植被 的 主要 气候 因素 。 成 方 妍 
等 . 引 分 析 气 象 因子 与 NDVI 相关 性 发 现 气象 因子 与 区 域内 不 同 土地 覆 被 类 型 NDVI 相关 性 不 显著 。 

近年 来 也 有 一 些 有 关 榆 林 及 周边 地 区 植被 指数 的 研究 , 李 登 科 等 "1 使 用 1981 一 2003 年 23a 长 序列 资料 
对 陕 北 长 城 沿线 风沙 区 植被 指数 利用 趋势 线 分 析 方 法 ,分 析 了 区 域 平 均 NDVI 的 变化 趋势 ,发 现 陕 北 长 城 沿 
线 风 涉 区 植被 覆盖 状况 尽管 有 波动 起 伏 ,但 整体 在 持续 转 好 ,年 均 NDVI 增加 了 10.62%。 并 利用 简单 相关 系 
数 的 方法 ,对 长 城 沿线 每 个 县 分 析 了 年 . 季 平 均 NDVI 与 年 份 的 相关 系数 和 和 斜率 。 刘 静 等 5 呈 利 用 NDVI 和 像 
元 分 解 模型 ,建立 了 毛乌素 沙漠 植被 覆盖 度 遥 感 定量 模型 ,并 对 1990 年 和 2007 年 2 个 时 期 植被 覆盖 度 进行 
了 等 级 划分 ,利用 转移 矩阵 的 方法 分 析 植 被 覆盖 度 的 变化 趋势 。 周 淑琴 轨 利 用 SPOT NDVI 数据 ,研究 了 不 同 
等 级 植被 间 的 变化 过 程 及 NDVI 空间 异 质 性 。 周 淑琴 等 5 应 用 空间 统计 学 和 经 典 统计 学 方法 ,研究 了 毛 乌 
素 沙漠 1998 一 2013 年 的 植被 空间 自 相 关 分 布 模式 和 发 展 趋势 及 气象 因素 的 影响 。 王 静 开 等 "利用 偏 相 关 
分 析 方 法 ,分 析 毛 乌 素 沙 漠 植 被 年 平均 物候 与 气温 和 降水 的 关系 。 这 些 研 究 采 用 的 植被 指数 的 时 间 序 列 、 空 
间 尺 度 .研究 方法 各 不 相同 。 郑 亚 云 55 利 用 2000 一 2014 年 MODIS 数据 ,采用 波段 重 运算 .最 大 化 合成 .时 间 
序列 分 析 线性 回归 分 析 方 法 研究 榆林 不 同 尺度 NDVI 的 变化 差异 ,同时 选取 空间 上 的 一 些 点 对 高 程 .坡度 及 
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气候 因子 与 NDVI 的 关系 做 定量 研究 。 刘 登 峰 等 ' 叶 基于 MODIS NDVI 遥感 数据 分 析 了 2000 一 2013 年 毛乌素 
沙漠 南部 植被 的 变化 状况 ,总 结 出 NDVI 年 均值 呈现 出 上 升 趋势 ,年 增长 率 为 2.69x10-?/a, 还 分 析 了 NDVI 逐 
年 的 级 别 转移 状况 。 

目前 对 NDVI 时 空 变化 分 析 , 有 些 是 基于 面 平 均 分 析 的 ,有 些 是 利用 一 元 线性 回归 趋势 线 方法 ,基于 像 元 
尺度 进行 分 析 。 目 前 对 NDVI 影响 因素 的 研究 涉及 以 下 3 个 方面 :对 研究 区 域 植被 指数 取 面 平均 值 与 降水 \ 
气温 等 气象 要 素 进行 简单 相关 分 析 或 回归 分 析 ; 在 像 元 尺度 的 简单 相关 分 析 ; 在 研究 区 域 选取 气象 站 点 周转 
的 NDVI 值 与 气象 因子 进行 简单 相关 分 析 和 多 元 回归 分 析 。 但 是 尚未 有 研究 对 某 一 特定 区 域 进 行 NDVI 像 元 
尺度 的 多 元 回归 分 析 , 并 定量 区 分 气象 因素 与 人 类 活动 对 植被 指数 的 影响 。 进 行 NDVI 像 元 尺度 的 分 析 , 可 
以 准确 分 析 研 究 区 域内 每 个 地 点 的 时 空 变化 及 气象 因素 对 小 范围 的 影响 ,空间 连续 性 和 异 质 性 能 够 更 好 的 体 
现 和 分 析 ; 利 用 多 元 回归 分 析 , 可 以 同时 考虑 多 个 气象 因素 对 植被 指数 的 影响 ,并 能 准确 分 析 气 象 因 素 对 植被 
指数 的 影响 及 贡献 ,实现 气象 因素 与 人 类 活动 对 植被 指数 影响 的 定量 区 分 ,比较 准确 的 基于 像 元 尺度 预测 未 
来 植被 覆盖 状况 。 所 以 ,本 研究 将 二 者 结合 起 来 ,对 陕西 榆林 植被 指数 进行 分 析 和 研究 。 

鉴于 像 元 尺度 研究 的 重要 性 ,本 文 以 陕西 省 榆林 市 为 研究 区 域 ,对 2000 一 2016 年 的 17 年 归 一 化 差分 植 
被 指数 在 像 元 尺度 进行 时 空 变化 分 析 ,然后 基于 像 元 尺度 利用 复 直 线 回 归 分 析 方法 ,对 植被 动态 变化 的 人 为 
因素 与 气象 因素 进行 定量 分 析 ,并 对 NDVI 进行 预测 。 


1 研究 区 域 概况 


榆林 市 位 于 陕西 省 最 北部 ,经 纬度 范围 是 107°%28'E 一 111°15 禄 ,36°57'N 王 39°35'N。 东 临 黄河 与 山西 省 
相 望 , 西 接 宁 夏 回 族 自 治 区 、 甘 肃 省 , 北 邻 内 蒙古 自治 区 , 南 接 陕 西 省 的 延安 市 。2015 年 底 , 全 市 常住 人 口 
340.11 万 人 。 地 域 东西 长 385km ,南北 宽 263km ,总 土地 面积 43578km*'”。 


107° 108° 109° 110° 111° 112°E 


图 1 榆林 市 位 置 和 地 形 图 
Fig.l Location and topographic map of Yulin City 


榆林 全 境地 势 呈 西北 高 、. 东 南 低 , 平 均 海拔 为 1300m ,地 貌 大 体 以 长 城 为 界 , 北 部 为 风沙 草滩 区 , 占 总 面积 
的 42% ,南部 为 黄土 丘陵 沟 蜜 区 , 占 总 面积 的 58% 。 榆 林 我 国 东部 季风 气候 和 西北 干旱 大 陆 性 气候 的 过 渡 地 
带 , 也 是 毛乌素 沙漠 南 缘 与 陕 北 黄土 高 原 的 过 渡 地 带 , 这 也 决定 了 榆林 生态 环境 的 脆弱 性 。 榆 林 属 暧 温带 半 
湿润 气候 向 半 干 旱 气 候 的 过 渡 区 ,属于 温带 大 陆 性 气候 ;从 东南 到 西北 随 着 距 海 距离 的 增加 ,植被 表现 出 明显 
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的 地 带 性 :温带 落叶 林 一 草原 一 鞠 澳 草原 一 攻 漠 的 过 渡 和 演 蔡 。 


2 数据 来 源 与 处 理 


本 文采 用 来 自 MODISZTerra 网 站 提供 的 NDVI 遥感 数据 ,数据 集 全 称 为 MODIS/Terra Vegetation Indices 
16-Day L3 Global 250m SIN Grid V005 ,简称 MOD13Q1。MOD13Q1 数据 是 16d 合成 ,空间 分 辨 率 为 230m ,数据 
文件 包含 了 12 个 字段 的 数据 ,NDVI 的 有 效 值 在 -1 到 1 之 间 !5 。 

本 研究 中 首先 对 下 载 的 影像 应 用 MRT( MODIS Reprojection Tool) 软件 进行 批量 裁剪 .投影 等 操作 ?由 于 本 
研究 只 对 8 月 份 NDVI 影像 进行 分 析 , MOD13Q1 数据 8 月 份 有 2 幅 影 像 , 得 到 后 的 影像 进行 最 大 化 合成 处 
理 。 这 样 每 年 得 到 一 幅 影像 ,从 2000 一 2016 年 共 17 幅 影像 。 


3 研究 方法 


3.1 一 元 线性 回归 趋势 线 分 析 

一 元 线性 回归 趋势 线 分 析 是 对 一 组 随时 间 而 变化 的 变量 进行 回归 分 析 的 方法 s 该 方法 能 够 分 析 研 究 区 
域 每 个 像 元 的 变化 趋势 ,以 此 来 模拟 植被 的 绿 度 变 化 率 ( Greenness Rate of Change ,CRC》, GRC 被 定义 为 某 时 
间 段 内 的 季节 合成 归 一 化 植被 指数 (Seasonally Integrated Normalized Difference Vegetation Index,SINDVT) 年 际 
变化 的 线性 回归 方程 的 斜率 :21 。 本 文通 过 每 个 像 元 17 年 的 NDVI, 可 以 模拟 该 像 元 的 NDVI 在 这 17 年 间 
的 变化 趋势 ,并 估计 变化 幅度 ,回归 直线 斜率 (slope) 采 用 最 小 二 乘法 求 得 ; 其 计算 公式 如 下 : 


17 x 3 x NDVT, — (3) (5 NDVE) 
slope = ee ~ Ps (1) 
17X XY - (之 
式 中 ;i 为 1 一 17 年 中 的 第 i 年 ; NDVI 为 第 i 年 的 NDVT 值 ; slope 为 趋势 线 的 斜率 。slope>0 时 ,说 明 NDVI 在 
这 17 年间 呈现 增加 的 趋势 ;slope=0 时 ,说 明 NDVI 没有 变化 ;slope<0 时 ,NDVI 呈现 减少 的 趋势 。 
3.2 ”相关 分 析 检 验 
自然 界 中 的 许多 现象 之 间 存 在 着 一 定 的 联系 ,它们 之 间 既 不 是 确定 的 函数 关系 ,也 不 是 完全 没有 关系 。 


相关 分 析 就 是 研究 两 个 或 多 个 随机 变量 之 间 的 联系 “”。 其 计算 公式 如 下 : 


> [ (x; -了 X) (y; -7)] 
一 (2 


> (x; -人 之 (=D 

式 中 : mn 为 样本 数 ; X 为 变量 x 的 均值 ; 了 为 变量 y 的 均值 ; r, 为 变量 x 和 变量 y 的 相关 系数 。 

研究 对 象 是 NDVI 序列 和 时 间 序 列 , 所 得 到 的 结果 为 每 一 像 元 对 应 的 NDVI 值 与 年 份 的 回归 系数 ,得 到 的 
相关 系数 态 " 可 以 用 来 检验 植被 生长 状况 的 变化 趋势 。 如 果 相 关系 数 为 正 值 ,表明 植被 覆盖 度 呈 现 增加 的 趋 
势 ; 反之 就 表示 该 区 域 植被 禾 盖 度 呈现 降低 的 趋势 。 如 果 r, 值 通过 了 0.05 的 显著 性 水 平 (P<0.05 ) 检验 ,说 
明 NDVI 具有 显著 的 变化 趋势 。 
3.3” 复 直线 回归 分 析 

本 文 是 进行 NDVI 与 温度 、 降 水 之 间 的 复 直 线 回 归 分 析 , 回 归 方 程 为 ; 

z=atbx+cey (3) 


参数 计算 公式 为 : 


b _ XY yy 。 了 (4) 
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TTT 0: 

Ee (5) 

a=z—bx—cy (6) 

式 中 ,x 、yY 、z 为 x 、y 、z 系 列 的 均值 ; go,、o，、o, 为 系列 的 均 方差 ;7 、r, 、r, 是 x* 和 y 、x 和 z、y 和 z 的 相 

NDVI 模拟 残 差 计算 公式 : 

residual = NDVI-NDVI, 全 

式 中 ,residual 为 NDVI 复 直 线 回 归 方程 残 差 ;NDVI 为 NDVI 时 间 序 列 数据 集 ; NDVI, 是 根据 复 直 线 回归 模型 
预测 的 NDVI 值 。 

4 结果 分 析 


4.1 ”NDVI 的 时 空 变化 特征 

通过 一 元 线性 回归 趋势 线 分 析 得 到 2000 一 2016 年 8 月 份 逐 像 元 NDVI 空间 分 布 坡度 图 (图 2) 和 不 同 级 
别 slope 面积 ( 表 1) ,图 2 中 绝 大 部 分 区 域 呈 现 绿 色 ,说 明 大 部 分 区 域 植被 指数 是 增加 的 ,增加 率 在 0 一 0.03/a 
占 多 数 ;有 些 区 域 达 到 了 0.04 以 上 ,最 大 值 达到 了 0.044 ;植被 指数 减 小 的 区 域 主 要 零星 分 布 在 榆林 市 的 西部 
及 西北 部 ,这 些 区 域 是 毛乌素 沙漠 的 南 缘 ,主要 涉及 到 定 边 上 县、 靖 边 县 \ 横 山区 和 榆 阳 区 。 最 小 值 为 -0.054/a， 
平均 值 为 0.0102/a ,标准 差 为 0.0058。 


40°N 
图 例 


slope 


-0.05 -0.03 -0.01 0 0.01 


39° 


38° 


37° 
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图 2 2000 一 2016 年 榆林 NDVI 一 元 线性 回归 坡度 分 布 图 
Fig.2 Distribution of NDVI one-dimensional linear regression slope in Yulin from 2000 to 2016 


表 1 为 榆林 植被 指数 一 元 线性 回归 分 析 变 化 表 , NDVI 减 小 的 面积 为 1537.44km?, 仅 占 榆林 总 面积 
3.57% ,其 中 slope 在 0 和 -0.1 之 间 的 区 域 为 1452.81km2 , 占 榆林 市 总 面积 的 3.37%;NDVI 增加 的 面积 ; 
41609.76km? , 占 榆林 总 面积 的 96.44% ,其 中 slope 在 0 一 0.02/a 之 间 的 区 域 占 93.63%。 说 明 近 些 年 来 榆林 
NDVI 呈现 增加 的 趋势 , 且 增 加 值 基本 都 在 0 与 0.02/a 之 间 ,榆林 植被 覆盖 度 逐 渐 增 加 。 

4.2 ”NDVI 变化 程度 评价 
为 了 定量 的 分 析 研 究 区 域 NDVI 的 变化 程度 ,对 逐年 的 NDVI 影 像 和 时 间 序 列 进行 相关 系数 分 析 ,得 到 每 
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一 像 元 的 NDVI 值 与 年 份 的 回归 系数 , 即 相关 系数 7 ,将 该 系数 进行 0.05 的 显著 性 水 平 检验 ,得 到 的 变化 趋 
势 分 布 图 见 图 3。 


表 1 2000 一 2016 年 榆林 NDVI 一 元 线性 回归 分 析 变 化 表 


Table 1 NDVI variance linear regression analysis of changes in the table in Yulin from 2000 to 2016 


NDVI 变化 趋势 /a 积 比例 NDVI 变化 趋势 /a 面积 比例 
NDVI variance tread Area/ km? Proportion/% NDVI variance tread Area/km? Proportion/ % 
slope<—0.05 0.32 0.00 0<slope<0.01 18349.61 42:53 
-0.05 <slope<-0.04 0.97 0.00 0.01<slope<0.02 22046.36 51.10 
-0.04<slope<-0.03 2.58 0.01 0.02<slope<0.03 1194.52 2.77 
-0.03<slope<-0.02 12.13 0.03 0.03 <slope<0.04 18.57 0.04 
-0.02 三 slope<-0.01 68.64 0.16 slope 宇 0.04 0.70 0.00 

-0.01<slope<0 1452.81 3.37 


40°N N 
和 
图 例 
39° 显著 性 检验 
国 国 显著 降低 
降低 但 不 显著 
增加 但 不 显著 
加 央 显著 增加 
38° 
Be 
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图 3 2000 一 2016 年 NDVI 动态 变化 显著 性 水 平 检 验 图 
Fig.3 The significance test of NDVI dynamic change from 2000 to 2016 


如 果 x; 值 大 于 0.05 的 显著 性 水 平 (P<0.05)roo =0.482 , 则 认为 NDVI 增加 或 碱 小 趋势 显著 。 图 3 显示 ， 
绝 大 多 数 区 域 通过 了 0.05 的 显著 性 水 平 检验 ,显著 增加 的 区 域 占 研 究 区 域 总 面积 的 80.72% ,增加 但 不 显著 
的 区 域 占 15.71%; 显 著 降 低 的 区 域 仅 占 0.43% ,零星 分 布 在 西北 部 毛乌素 沙漠 处 ,降低 但 不 显著 的 区 域 占 3. 
149%, 零 星 分 布 在 西部 的 定 边 县 和 靖 边 县 。 所 以 ,榆林 的 NDVI 近 些 年 来 呈现 显著 增加 的 趋势 , 呈 显 著 增 加 趋 
势 的 面积 占 80.72% 。 
4.3 ”2000 年 与 2016 年 比较 分 析 

将 榆林 2016 年 8 月份 NDVI 影像 减 去 2000 年 8 月 份 NDVI 影 像 ,得 到 两 个 年 份 的 差 值 分 布 图 (图 4) ,并 
统计 出 不 同 差 值 范围 所 占 的 面积 ( 表 2) 。2016 年 与 2000 年 相 比 ,大 部 分 区 域 NDVI 增 加 值 在 0 一 0.4 之 间 。 

2016 年 NDVI 值 比 2000 年 小 的 区 域 仅 占 研究 区 域 的 3.14% ,为 1354.12km ;其 中 小 于 0.2 的 仅 有 42.45 
km , 占 总 面积 的 0.10% ;在 0.2 与 0 之 间 的 为 1311.67km?, 占 榆林 总 面积 的 3.04%。2016 年 NDVI 值 比 2000 
年 增加 的 区 域 占 研 究 区 域 的 96.86% ,达到 了 41793.09km? ; 其 中 在 0 与 0.2 之 间 的 区 域 是 19841.28km2 , 占 研 
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究 区 域 总 面积 的 45.99% ;在 0.2 与 0.4 之 间 的 有 21505.8$km2 , 占 总 面积 的 49.84% ;在 0.4 和 0.6 之 间 的 为 
438.71km” , 仅 占 1.02% ;大 于 0.6 的 区 域 仅 有 7.25km” , 仅 占 总 面积 的 0.02%。 说 明 2016 年 与 2000 年 相 比 , 榆 
林 市 绝 大 部 分 区 域 植被 指数 都 增加 了 , 且 增 加 幅度 基本 在 0 一 0.4 之 间 。 


40°N 上 图 例 
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国 国 -0.76 一 -0.6 
国 一 -0.6 一 -0.4 
-0.4—0.2 
-0.2 一 0 
0 一 0.2 
39° 0.2 一 0.4 
0.4—0.6 
0.6—0.78 
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图 4 2016 年 与 2000 年 8 月 份 NDVI 差 值 分 布 图 
Fig.4 NDVI difference distribution in August between 2016 and 2000 


表 2 2016 年 与 2000 年 榆林 8 月 份 NDVI 差 值 分 布 表 
Table 2 NDVI difference distribution table in in August between 2016 and 2000 


NDVI 差 值 面积 比例 NDVI 差 值 面积 比例 
i NDVI difference Area/ km? Proportion/ % NDVI difference Area/ km? Proportion/ % 
= NDVI<-0.6 1.45 0.00 0< NDVI<0.2 19841.28 45.99 
-0.6<NDVI<-0.4 4.13 0.01 0.2< NDVI<0.4 21505.85 49.84 
-0.4<NDVI<-0.2 36.87 0.09 0.4< NDVI<0.6 438.71 1.02 
-0.2<NDVI<0 1311.67 3.04 NDVI=0.6 7:25 0.02 


为 了 分 析 2000 一 2016 年 不 同 级 别 植被 指数 的 变化 情况 ,参照 周 淑琴 中 的 研究 ,将 榆林 植被 序列 分 为 6 个 
等 级 :1 级 CNDVI<0.2) 无 植被 区 2 级 (0.2<NDVI<0.3) 极 低 植 被 履 盖 区 、3 级 (0.3<NDVI<0.4) 低 植被 覆盖 
区 .4 级 (0:4<NDVI<0.5) 中 等 植被 覆盖 区 5 级 (0.5<NDVI<0.6) 高 植被 覆 羡 区 .6 级 (NDVI>0.6) 密集 植 被 覆 
盖 区 。 做 出 2000 年 植被 指数 等 级 向 2016 年 植被 指数 等 级 的 转移 矩阵 ( 表 3) 。 

2000 一 2016 年 榆林 市 植被 等 级 面积 转移 矩阵 表 ( 表 3) 中 ,对 角 线 上 的 数据 是 未 发 生 转 变 的 植被 ,对 角 线 
右上 方 是 正 向 转变 的 植被 等 级 面积 ,左下 方 是 逆向 转变 的 植被 等 级 面积 。 显 然 右上 方 的 面积 明显 大 于 左下 方 
的 面积 ,这 和 作 图 分 析 得 到 的 结果 相符 合 。 

4.4 植被 动态 变化 气象 因素 与 人 为 因素 的 定量 分 析 

为 了 逐 像 元 定量 分 析 气 象 因素 和 人 类 活动 对 植被 的 影响 ,另外 由 于 气象 因子 存在 空间 不 均匀 性 ,选取 榆 
林 区 域内 及 其 周边 的 13 个 气象 站 点 资料 ,采用 空间 插值 的 方法 来 获得 每 一 个 像 元 的 降水 和 气温 数据 。 选 用 
了 榆林 、 定 边 . 靖 边 .横山 ` 绥 德 .神木 . 吴 旗 、 延 安 . 兴 县 、 离 石 .河曲 .盐池 、 东 胜 共 13 个 站 点 2000 一 2015 年 的 
降水 和 和 气温 资料 。 由 于 植被 指数 与 温度 .降水 之 间 存 在 一 定 的 滞后 关系 ,所 以 为 了 准确 分 析 NDVI 与 温度 和 
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降水 之 间 的 定量 关系 ,必须 考虑 这 种 滞后 效应 。 在 选取 月 降水 和 气温 数据 时 ,利用 泰 森 多 边 形 方法 将 点 降水 
和 气温 分 配 到 面 上 ,分别 求 出 榆林 NDVI 与 6 月 7 月 8 月 .6 月 和 7 月 平均 值 7 月 和 8 月 平均 值 6 月 7 月 8 
月 这 3 个 月 平均 值 的 相关 系数 ,并 进行 了 1 分 布 检验 ,发 现 榆 林 NDVI 与 7 月 份 平均 气温 和 7 月 份 总 降水 量 t 
值 最 大 , 均 通 过 了 a =0.05 的 置信 和 度 检验 。 所 以 本 文选 择 了 2000 一 2015 年 13 个 站 点 7 月 份 平均 气温 和 7 月 
份 总 降水 量 作 为 气温 和 降水 的 研究 序列 。 


表 3 2000 一 2016 年 榆林 市 植被 等 级 面积 转移 矩阵 表 / km? 
Table 3 Vegetation grade area transfer matrix table in Yulin from 2000 to 2016 


2016 年 

等 级 1 级 2 级 3 级 4 级 5 级 6 级 总 计 

Grades Level 1 Level 2 Level 3 Level 4 Level 5 Level-6 Total 
2000 年 1 级 155.47 1677.96 1344.72 726.12 365.87 80.89 4351.03 
2 级 52.28 1278.88 3511.10 3833.10 5121.33 1905.58 15702.27 
3 级 9.82 140.68 1216.95 2980.22 5800.50 $5830.18 15978.36 
4 级 1.99 20.26 178.19 738.60 1615.96 2933.65 5488.65 
5 级 1,77 5.98 20.62 70.70 273.75 781.89 1154.72 
6 级 2:52 3.43 6.80 10.24 38:47 406.63 468.10 
总 计 Total 223.85 3127.19 6278.39 8358.97 13215.88 11938.83 43143.12 


把 NDVI 的 2000 一 2015 年 时 间 序 列 和 所 采用 的 气温 和 降水 数据 进行 复 直 线 回归 分 析 , 求 出 a、b ce 参数 ， 
再 计算 出 NDVI 与 温度 .降水 之 间 的 回归 方程 ,最 后 利用 温度 和 降水 序列 得 出 NDVI 的 预测 值 。 该 预测 值 是 气 
候 因 素 对 NDVI 的 贡献 部 分 ,用 原始 的 NDVI 减 去 NDVI 的 预测 值 即 可 得 到 残 差 ,该 残 差 即 为 人 为 因素 对 
NDVI 的 影响 和 其 他 不 确定 因素 的 影响 。 

图 5 是 将 2015 年 榆林 面 气温 和 降水 数据 进行 复 直 线 回 归 分 析 得 到 的 NDVI 值 ,该 图 东部 和 东南 部 区 域 
NDVI 多 在 0.4 以 上 ,西部 和 西北 部 的 风沙 草原 区 NDVI 值 较 小 ,说 明 东 部 和 东南 部 的 森林 草原 区 自然 条 件 
(地 形 地貌 .气候 水文、 土壤 ) 更 适合 于 植被 的 生长 ;西部 和 西北 部 由 于 位 于 毛乌素 沙漠 的 边缘 ,自然 条 件 恶 
劣 ,植被 自我 更 新 和 自我 生长 能 力 较 差 , 表 4 是 与 图 5 相对 应 的 不 同 NDVI 贡献 值 的 面积 及 比例 。99% 以 上 
的 区 域 NDVI 为 正 值 ,其 中 0.3<NDVI<0:4 的 面积 为 9220.72km? , 占 到 总 面积 的 21.46% ;0.4 生 NDVI<0.5 的 面 
积 为 18612.14km” , 占 总 面积 的 43.32% ;0.5 和 NDVI<0.6 的 面积 为 8763.92km, 占 总 面积 的 20.40%。 三 者 总 
和 占 到 了 总 面积 的 85% 以 上 ,说 明 在 没有 人 类 活动 影响 下 , 现 阶 段 榆 林 的 NDVI 一 般 在 0.3 一 0.6 之 间 。 而 
2015 年 的 实际 观测 值 0.3 志 NDVI<0.4 占 到 总 面积 的 26.89% ;0.4 和 NDVI<0.5 占 到 总 面积 的 31.10% ;0.5 三 
NDVI<0.6 占 到 总 面积 的 .22.90% ;0.2 入 NDVI<0.3 的 比例 为 12.22%;0.6 入 NDVI<0.7 的 比例 为 5.32%。 实 际 
观测 值 与 预测 值 的 主要 差别 为 介 于 0.4 和 0.5 之 间 的 区 域 面积 减 小 了 , 减 小 了 12.22% ,而 0.6 和 0.7 之 间 的 区 
域 增 加 了 5 增加 了 3.60% 。 同 时 还 可 以 看 出 人 类 活动 使 得 植被 指数 变 得 更 加 均匀 ,中 等 植被 覆盖 区 面积 减少 
了 ,而 密集 植被 覆盖 区 面积 增加 了 ,这 是 近年 来 榆林 市 开展 封山育林 、 退 耕 还 林 、 退 牧 还 草 等 各 项 植被 恢复 和 
保护 措施 的 结 


表 4 气象 因素 (气温 和 降水 ) 对 榆林 市 2015 年 8 月 NDVI 贡献 值 分 布 表 
Table 4 Meteorological factors (temperature and precipitation) to NDVI contribution distribution table for Yulin City in August 2015 


NDVI 贡献 值 面积 比例 NDVI 贡献 值 面积 比例 
NDVI contribution Area/ km? Proportion/ % NDVI contribution Area/ km’ Proportion/ % 
NDVI<0 21.73 0.05 0.4< NDVI<0.5 18612.14 43.32 
0<NDVI<0.1 27.37 0.06 0.5< NDVI<0.6 8763.92 20.40 
0.1< NDVI<0.2 449.98 1.05 0.6< NDVI<0.7 737.84 1.72 
0.2< NDVI<0.3 5030.20 11.71 0.7< NDVI<0.8 104.75 0.23 
0.3<NDVI<0.4 9220.72 21.46 NDVI=0.8 0.06 0.00 
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图 5 气象 因素 (气温 和 降水 ) 对 榆林 市 2015 年 8 月 NDVI 贡献 分 布 图 
Fig.5 Meteorological factors (temperature and precipitation) to NDVI contribution distribution for Yulin City in August 2015 


为 了 定量 化 分 析 人 类 活动 对 植被 的 影响 ,根据 复 直 线 回归 分 析 , 用 原始 观测 的 NDVI 减 去 NDVI 的 预测 
值 , 即 计 算出 了 残 差 , 残 差 可 以 认为 是 人 为 因素 对 NDVT 所 贡献 的 部 分 (图 6 和 表 5)。 人 类 活动 对 植被 指数 
有 消减 的 区 域 占 榆林 总 面积 的 45.04% ,主要 分 布 在 榆林 最 北部 的 府 谷 县 .榆林 南部 .榆林 西部 等 区 域 ( 图 6)。 
人 类 活动 对 植被 指数 有 提高 的 区 域 占 榆 林 总 面积 的 54:96% ,但 有 46.64% 的 区 域 的 贡献 值 介 于 0 和 0.1 之 间 ， 
说 明 2015 年 这 些 区 域 的 人 类 活动 使 得 植被 指数 增加 了 ,植被 覆盖 度 增加 。2015 年 榆林 市 统计 年 鉴 显示 ， 
在 植被 指数 增加 的 区 域 森林 履 盖 率 ` 第 一 产业 增加 值 .年末 常用 耕地 面积 比例 、 水 保 治理 面积 .旱涝保收 面积 
都 明显 大 于 植被 指数 降低 的 地 区 。 榆林 市 从 1999 年 开始 积极 响应 国家 号 召 ,逐渐 开展 退耕 还 林 荒山 荒地 造 
林 、 封 山 育 林 等 工程 。 这 些 退 耕 还 林 还 草 的 生态 工程 建设 促进 了 植被 的 恢复 ,但 是 这 些 措 施 在 不 同 地 区 的 实 
施 规模 和 实施 效果 是 存在 差异 的 而 NDVI 的 变化 直接 反映 了 植被 恢复 的 效果 ,所 以 此 方法 可 以 作为 监察 封 
山 育 林 、 退 耕 还 林 < 退 牧 还 草 的 政策 落实 和 实施 效果 的 参考 ,鉴于 统计 回归 分 析 结 果 存 在 的 误差 ,在 应 用 中 应 
予以 注意 。 男 外 在 现 阶 段 人 类 活动 变化 幅度 不 大 的 情况 下 ,可 以 根据 未 来 的 气温 和 降水 预测 未 来 植被 状况 的 
理论 值 。 


表 5 人 类 活动 对 榆林 2015 年 8 月 NDVI 的 贡献 值 分 布 表 
Table 5 Human activity to NDVI contribution distribution table for Yulin City in August 2015 


NDVI 贡献 值 面积 比例 NDVI 贡献 值 面积 比例 
NDVI contribution Area/ km? Proportion/ % NDVI contribution Area/ km? Proportion/ % 
NDVI<0 19349.17 45.04 0.3< NDVI<0.4 56.99 0.13 
0<NDVI<0.1 20036.08 46.64 0.4< NDVI<0.5 17.07 0.04 
0.1< NDVI<0.2 3285.51 7.65 NDVI=0.5 1:72 0.00 
0.2< NDVI<0.3 210.74 0.49 
5 结论 


本 文 利用 MODIS/Terra NDVI 时 间 序 列 数 据 对 2000 一 2016 年 陕西 榆林 生长 季 (8 月 份 ) 植 被 指数 进行 时 
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图 6 人 类 活动 对 榆林 2015 年 8 月 NDVI 贡献 分 布 图 
Fig.6 Human activity to NDVI contribution distribution for Yulin City in August 2015 


空 变化 特征 及 影响 因素 分 析 , 主 要 得 出 以 下 结论 : 

(1) 对 2000 一 2016 年 榆林 8 月 份 NDVI 影像 进行 时 空 变化 分 析 表 明 ,榆林 有 96.44% 的 区 域 植被 指数 是 
增加 的 ,增加 率 在 0 一 0.02/a 之 间 的 区 域 占 榆 林 的 93:63%, 显著 降低 的 区 域 仅 占 0.43% ,零星 分 布 在 西北 部 
毛乌素 沙漠 处 , 呈 显 著 增 加 趋势 的 面积 占 到 80.72% 。2016 年 与 2000 年 相 比 ,榆林 绝 大 部 分 区 域 植被 指数 都 
增加 了 , 且 增 加 幅度 基本 在 0 一 0.4 之 间 。 

(2) 为 了 逐 像 元 定量 分 析 气 象 因素 和 类 活动 对 植被 的 影响 , 复 直 线 回 归 分 析 显 示 , 气 象 因 素 对 植被 生 
长 演化 起 到 了 促进 作用 ,人 类 活动 使 得 植被 指数 等 级 变 得 更 加 均匀 。 气 象 因素 对 植被 指数 的 贡献 东部 和 东南 
部 区 域 的 森林 化 草原 区 NDVI 多 在 0.4 以 上 ,西部 和 西北 部 的 风沙 草原 区 NDVI 的 贡献 较 小 ,说 明 东 部 和 东南 
部 区 域 自然 条 件 ( 地 形 .地貌 . 气 候 、 水 文 .土壤 ) 更 适合 于 植被 的 生长 ,西部 和 西北 部 由 于 位 于 毛乌素 沙漠 的 
边缘 ,自然 条 件 恶劣 ?植被 自我 更 新 和 自我 生长 能 力 较 差 。 

(3) 人 类 活动 表现 为 抑制 植被 生长 演化 的 区 域 占 榆林 总 面积 的 45.04% , 主要 分 布 在 榆林 市 最 北部 的 府 谷 
县 .榆林 南部 < 榆林 西部 等 区 域 。 人 类 活动 对 植被 指数 有 提高 的 区 域 占 榆 林 总 面积 的 54.96% , 说明 有 一 半 多 
以 上 的 区 域 ,人 类 活动 对 植被 生长 起 到 了 促进 作用 ,这 些 地 区 的 封山育林 、 退 耕 还 林 、 退 牧 还 草 等 措施 的 实施 
效果 较 好 。 在 现 阶段 人 类 活动 变化 幅度 不 大 的 情况 下 ,可 以 根据 未 来 的 气温 和 降水 预测 未 来 植被 状况 的 理论 
值 , 为 区 域 生 态 修 复 和 环境 保护 提供 参考 。 
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